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ABSTRACT 

We investigate the occurrence of crystalline silicates in oxygen-rich evolved stars across 
a range of metallicities and mass-loss rates. It has been suggested that the crystalline 
silicate feature strength increases with increasing mass-loss rate, implying a correlation 
between lattice structure and wind density. To test this, we analyse Spitzer IRS and 
Infrared Space Observatory SWS spectra of 217 oxygen-rich asymptotic giant branch 
and 98 red supergiants in the Milky Way, the Large and Small Magellanic Clouds and 
Galactic globular clusters. These encompass a range of spectral morphologies from the 
spectrally-rich which exhibit a wealth of crystalline and amorphous silicate features 
to 'naked' (dust-free) stars. We combine spectroscopic and photometric observations 
with the GRAMS grid of radiative transfer models to derive (dust) mass-loss rates and 
temperature. We then measure the strength of the crystalline silicate bands at 23, 28 
and 33 /im. We detect crystalline silicates in stars with dust mass-loss rates which 
span over 3 dex, down to rates of ~10 -9 M yr _1 . Detections of crystalline silicates 
are more prevalent in higher mass-loss rate objects, though the highest mass-loss rate 
objects do not show the 23-/im feature, possibly due to the low temperature of the 
forsterite grains or it may indicate that the 23- fim band is going into absorption due 
to high column density. Furthermore, we detect a change in the crystalline silicate 
mineralogy with metallicity, with enstatitc seen increasingly at low metallicity. 

Key words: stars: AGB, post-AGB - circumstellar matter - dust, extinction - 
Magellanic Clouds - infrared: stars - radiative transfer 



1 INTRODUCTION 

Cool evolved stars in the later stages of active nuclear 
burning expel a large fraction of their mass through 
a slow, dense stellar wind at rates of 10~ 7 to 10~ 4 
M yr- ^Habing fc Olofssonll2003h . In low- to intermediate- 
mass stars (0.8-8 Mq) this occurs on the asymptotic giant 
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branch (AGB), while higher-mass stars (8-25 Mq) undergo 
this intense mass loss as red supergiants (RSG). In AGB 
stars it is believed that radial pulsatio ns levitate gaseous 
mate rial above the surface of the star l|Habing fc Olofssonl 
2003), resulting in a cooling outflow from which molec- 
ular species form and (at temperatures b elow ~1400 K) 
dust grains condense (|Hoefner et al.| [l998). A dust-driven 
wind develop s, propelled by radiation pressure from the stel- 
lar photons i|Gehrz fc WoolJ Il97ll : H owcn fc Willson|[l99ll : 
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iNorris et al.ll2012h . For RSGs, a slightly different me chanism 
may drive the dense outflow ([josselin fc Pledl2007l ). 

Chemistry in AGB atmospheres is affected by mixing 
of the stellar atmosphere and interior, known as dredge 
up. This process brings nuclear-fusion products, particularly 
carbon, to the stellar surface, enhancing their abundances. 
Chemistry in the cooling wind is dominated by the relative 
abundances of carbon and oxygen, as these combine at high 
temperatures to form carbon monoxide. The remaining free 
carbon or oxygen then determines whether the star produces 
C-rich dust (such as amorphous carbon and SiC) or O-rich 
dust (such as silicates and metal oxides). As we are inter- 
ested in silicate dust production, we focus this paper only 
on O-rich stars. 

Infrared (IR) spectroscopic observations of O-rich 
evolved stars show a plethora o f fea t ures due to 
molecules and dust JOmont et al.l 19931; Waters et al.l 



Table 1. LMC and SMC Parameters 



ll99rj ; ISvlvester et al.|[l999l ; iMolster et alj|2002al ) . Molecules, 
mainly CO, OH, H2O, and SiO, dominate the spectra at A < 
Wfim. At longer wavelengths, dust features become more 
important. Silicate dust dominates the IR spectra of the ma- 
jority of dust-producing O-rich AGB stars, and can be found 
in either amorphous or crystalline form. Amorphous silicates 
are most abundant, characterised by their broad, smooth 
features at 10 and 18 /im. Crystalline silicates exhibit sharp 
resonance features at 10 /im and longer wavelengths. The 
positions and shapes of these are very sensitive to composi- 
tional changes, lattice structure, and grain size and morphol- 
ogy, providing a mechanism to identify the chemical com- 



positio n and mineralogy of t he dust grai ns ( Molstc r et al 



2002bl; Ichiharaet al.l [20021: I Koike et all 120031 ; iMin et al 



20031 : IMolster fc Kemperll2005l ). Comparisons with labora- 



tory data have primarily identified these grains as Mg-rich 
olivines and pyroxenes (forsterite: Mg2 SiQ4 and enstatite : 
MgSiOa) with little or no iron content |Molster et al. 1999: 
Ide Vries et"aill201ll 1. 

The formation of crystalline silicates in the wind of 
evolved stars is not well understood: crystalline structures 
are the energetically more favourable atomic arrangement 
for silicates yet, in most cases, crystalline silicate dust is only 
a minor component of the circumstellar dust shell. The tran- 
sition from amorphous to crystalline silicate grains is a ther- 
ma l process, requiring temperatures in the region of 1040 
K (iHallenbeck et all ll99St (Fabian et all 120001 : [Speck et all 
l201lf ). However, there is little consensus on how silicates in 
stellar outflows gain sufficient energy for this transition to 
occur. In order to determine the physical conditions under 
which crystalline silicates form we require a better under- 
standing of the conditions surrounding the star (e.g. wind 
densities). The formation process of crystalline grains can be 
constrained by correlating the crystalline fraction with the 
dust or gas column density. A correlation with dust density 
suggests that annealing of amorphous silicate grains heated 
by radiation is probabl y the dominant means b y which crys- 
tals are manufactured ( Sogawa fc Kozasall 19991 ). while a cor- 
relation with gas density would suggest that direct conden- 
sation of crystalline silicates in the w ind will dominate (e.g. 
lTielendll998l : lGail fc Sedlmavrlll999l ). 

In Infrared Space Observatory (ISO) spectra of Milky 
Way (MW) giant stars, the spectral features due to the crys- 
talline silicates forsterite (Mg2Si04) and enstatite (MgSiOa) 
typically only appear around evolved stars if their mass- 



Parameter 


LMC 


Ref. 


SMC 


Ref. 


Distance, d (kpc) 


51 ±2 


1,2 


61 ± 2 


1,2 


Metallicity, Z (Zq) 


0.5 ±0.17 


3,4 


0.2 ±0.06 


3,4 


Inclination angle (°) 


34.7 ±6.2 


5 


68 ± 2 


6 


Gas-to-dust ratio ( <?) 


200 




500 




E(B-V) (mag) 


0.13 


7 


0.04 


8 



Refer enc es: (1) ICioni et al] ll2000h: (2) ISzewczvk et all 
J2009h; (3) iLuck et al.l dl998h ; (4) iMeixner etaL |201C ); (5) 
van der Marel fc CionH iboOli); f6) iGroenewegen feooof) ; (7) 
Massev et al] dl995h ; ~8l lHarris fc Zaritskvl ]2004l) . 



loss r ate is higher than a threshold value of ~10 5 Mq 
yr" 1 (|Cami et al.lll998l ; ISvivester et al.lll999T ). For example, 
the high-density winds of heavily enshrouded OH/IR stars 
have measured cr ystalline fractions of up to ~20 per cent o f 
the silicate mass jKemper et alj|200ll : ld"e Vries et alj|20ich . 
Conversely, lower mass-loss rate AGB stars, such as Miras, 
normally lack crystalline silicate dust features in their spec- 
tra. This suggests a correlation be t ween lattice order and 
wind density (e.g. ICami et al.lll998l : ISpeck et al.ll2008l ). Al- 
ternatively, contrast effects between amorphous and crys- 
talline silicate grains at different temperatures can mask the 
characteristic spectral features of the crystalline material 
(up to the 40 per cent mass fraction l evel) in infrared spec - 
tra of low mass- loss rate AGB stars (|Kemper et al.l [200lh . 
If the crystalline and amorphous silicates are in thermal 
contact with each other, contrast improves, and the detec- 
tion of smaller amounts (<40 per cent) of crystalline sili- 
cates becomes possible. This could explain observations of 
crystalline silicates in lower mass-loss rate objects, such as 
the recent observations of low-metallicity, low mass-loss rate 
evolved stars in Galactic globular clusters that show crys- 
talline silicates (ISloan et al] l20ld ; iMcDonald et al.l l2011al : 
iLebzelter et al]|2006l V ~ 

Observations of evolved stars in the metal-poor envi- 
ronments of the Magellanic Clouds with the Spitzer Space 
Telescope provide an ideal opportunity to explore the occur- 
rence of crystallinity and investigate how the O-rich dust 
condensation depends on the physical and chemical condi- 
tions of the envelope. Within the Large Magellanic Cloud 
(LMC) and Small Magellanic Cloud (SMC), we have a set 
of stars with similar characteristics: for instance, they are 
all at approximately the same distance from the Sun and 
are assumed to have a single metallicity within each galaxy. 
This allows us to measure luminosities and dust mass-loss 
rates from the observed spectral energy distributions. The 
parameters adopted in our calculations for the Magellanic 
Cloud stars are listed in Table Q] The dust-to-gas mass ra- 
tios of stars in the Magellanic Clouds are assumed to be 
lower than those of stars in the Solar nei ghbourhood and 
have a linear dependence with meta llicity ( van Loonll2000l : 
iMarshall et al.ll2004l ; Ivan Loonll2006l) . breaking the degener- 
acy between the dust column density and the gas density in 
the outflows of AGB stars. By carefully studying the depen- 
dence of crystallinity on the dust and gas mass-loss rate, we 
address the influences of dust density and gas density on the 
formation of crystalline grains. 

This paper is organised as follows: Section [2] describes 
the samples. In Section [3] we fit the spectral energy distri- 
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bution (SED) of each source and measure the strength of 
crystalline silicate emission bands. Average feature profiles 
are shown in Section[4] and the results are presented. Finally 
in Section[5]we discuss the implications for the onset of crys- 
tallinity and the change in mineralogy with metallicity. The 
conclusions are summarised in Section [S] 



2 THE SAMPLE 

Our sample contains 69 oxygen-rich AGB stars (O-AGB) 
and 76 RSG stars in the Magellanic Clouds which were 
observed spectroscopically with Spitzer, and 131 Galactic 
field O-AGBs and RSGs observed with either Spitzer or 
ISO. The Spitzer spectra cover a wavelength range of 5.2- 
37.2 fim while ISO spectra cover the 2.38-45.2 /im part of 
the spectrum. We combine this sample with 39 spectra from 
14 Galactic globular clusters to extend the low end of the 
metallicity range. 



2.1 LMC sample 

The LMC sample comprises 54 O-AGB and 6 RSG stars 
obser ved using the Infrared Spectrograph (IRS. lHouck et all 
I2004T ) on board Spitzer, either as p art of the SAGE-Spec 
legacy program jKemper et alj20ld ) , a spectroscopic follow- 
up to the Surveying the Agents of a Ga laxy's Evolution 
project (SAGE-LMC. iMeixner et alj|2006h . or drawn from 
the homogeneously-reduced archival IRS staring mode spec- 
troscopy within the SAGE-LMC footprint which are incor- 
porated into the SAGE-Spec data delivery. For a full de- 
scription of the original target selection, observing strategy 
and the techniques used in the data reduction for the SAGE- 
Spec l egacy program the reader is referred to lKemper et all 
<|2010h . 

For all sources in the Magellanic Clouds the asso- 
ciated broad-band photometry, including optical UBVI 
photome try from the Mag ellanic Clouds Photometric 
Survey dZaritsky et all l2004h . 2MASS JHK S photometry 
jSkrutskie et allhood ). mid-IR photometry (IRAC 3.6, 4.5, 
5.8, 8.0 and MIPS 24 /nm) and far-IR photometry (MIPS 
70 and 160 M m )i wa s compiled from the SAGE catalogue 
l|Meixner et al.ll2006l ). 

The O-AGB and RSG sources in the sample were se- 
lected b ased on the de c ision- tree classification scheme devel- 
oped bv lWoods et alj i|201ll ). Within the SAGE-Spec cata- 
logue, 197 objects have been classified according to their 
Spitzer IRS spectrum and SED with associated U, B, V, I, 
J, H, K s , IRAC and MIPS photometry. From this we found 
40 O-AGB stars, and 19 RSGs. To supplement this LMC 
sample, we applied th e decis ion-tree classification method- 
ology of IWoods et alj l|201ll ) to all the archival IRS ob- 
servations of point sources within the SAGE-LMC foot- 
print (Woods et al., in prep.), resulting in a further 148 
sources. As this work is primarily concerned with crystalline 
silicates, which become prominent at longer wavelengths 
(A > 20 /im), for inclusion in the sample we require that the 
sources were observed with the Long-Low (14.5-37.2 /xm) 
module on IRS. We limit our study to those spectra which 
were visually deemed to have a sufficient signal-to-noise for 
the identification of crystalline silicate features. Table[2]lists 
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Figure 1. Comparison of the bolometric magnitudes in the LMC, 
SMC and globular cluster samples. The 'classical' luminosity limit 
for AGB stars is indicated by the dotted line. Sources to the left 
of the line are AGB stars (blue) and to the right RSGs (red). 



the sources in the LM C sample, along with SAGE-Spec iden- 
tification (SSID) (see lKemper et alfoich . coordinates, As- 
tronomical Observation Request (AOR), bolometric magni- 
tudes (Mboi), periods and mass- loss rates determined from 
SED fitting with the GRAMS model grid (see Section |3~2|) . 
Only the first entries are shown for guidance, a full version 
is available in the electronic edition. 



2.2 SMC sample 

The SMC sample of 15 O-AGB and 16 RSG stars was com- 
piled f rom IRS staring mod e observations within the SAGE- 
SMC ((Gordon et al.l l201ll ) spatial coverage. Again, these 
were supplemented with optical, 2MASS, IR AC and MIPS 
photo metry from the SAGE-SMC ca talogue llGordon et al.1 
120111 ). and the S 3 MC catalogue (jBolatto et al.1 I2007T I 
Evolved stars were identified from th e SAGE-SMC mid- 
IR classifications of iBover et al.1 l|201ll ). with 57 candidate 
evolved stars found. IR photometry alone is insufficient to 
accurately separate the dusty O-AGB and RSG stars from 
contamination by young stellar objects (YSOs) with strong 
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Table 2. Details of the sources in the LMC sample for the first few entries. A full version is available in the electronic edition. For each 
source the SAGE-Spec identification (SSID), name, coordinates, Astronomical Observation Request (AOR), classification, bolometric 
magnitudes (M^oi), period and dust mass-loss rate (Mdust) is given. 



SSID 


Name 


RA 


Dec 


AOR 


Class 


M bo l a 


Period 


Mdust 






(J2000) 


(J2000) 








(d) 


(M @ yr- 1 ) 


4007 


IRAS 04407-7000 


04 40 28.49 


-69 55 12.70 


6077440 


O-AGB 


-6.68 


1179 6 


3.61E-08 


4 


SSTISAGE J044718.64-694220.5 


04 47 18.63 


-69 42 20.53 


22400256 


RSG 


-6.81 




5.01E-09 


4022 


MSX LMC 1132 


04 49 22.46 


-69 24 34.90 


10958336 


RSG 


-7.97 


361 d 


1.19E-08 


6 


SSTISAGE J044934.30-690549.2 


04 49 34.38 


-69 05 49.17 


24318720 


O-AGB 


-5.54 




2.59E-09 


4024 


IRAS 04498-6842 


04 49 41.38 


-68 37 51.00 


6076672 


RSG 


-7.86 


1280 d 


5.43E-08 


4031 


IRAS 04509-6922 


04 50 40.43 


-69 17 31.70 


6022400 


RSG 


-7.51 


1292 d 


2.85E-09 


8 


SSTISAGE J045128.56-695550.1 


04 51 28.58 


-69 55 49.90 


22400768 


O-AGB 


-3.82 


884 6 


2.09E-08 


4038 


IRAS 04516-6902 


04 51 28.94 


-68 57 49.50 


6020096 


O-AGB 


-5.48 


1090 d 


5.12E-09 


4054 


2MASS J04524318-7047371 


04 52 43.18 


-70 47 37.80 


33283584 


O-AGB 


-4.22 




7.90E-11 



(a) Bolometric m agnitudes for SS I D 1-197 are tak en from I Woo ds et al 1 toll) , all other values arc calculated as described in Sec. 13.1 
References: (b) ISoszvhski et alj [|2008l . l2009l) ; (c) lFraser et alj d2005l . l200Sh ; fd) iGroenewegen et all d2009h ; (e) iPoimanskH d2002Tl . 



oxygen-rich dust features, in particular those with a 10-[im 
silicate absorption feature. As with the LMC sample, we 
classify each object according to its Spitzer IRS spectrum 
based on its spectral features, SED, calculated bolometric 
luminosity, variability information and additional informa- 
tion from the literature, returning 27 O-AGB stars and 18 
RSGs. (A full spectral classification for the SMC is ongoing 
and will be presented in Ruffle et al in prep.). Again, sources 
observed solely with the IRS Short-Low setting or with a 
poor signal-to-noise ratio were excluded from the sample. A 
summary of the SMC sample can be found in Table [3] 

2.3 Galactic field sample 

To provide a solar metallicity comparison, our sample 
also includes ISO Sho rt Wavelength Spectrometer (SWS; 
Ide Graauw et all I1996T ) spectra of 79 O-AGB, 9 probable 
O-AGBs and 22 RSG stars in the Milky Way, covering a 
wavelength range of 2.38 - 45.2 /im. This Galactic field 
sample was compiled from the ISO/SWS databas^E which 
contains 1271 observations in full-scan mode. This database 
contains observations obtained for a wide range of individual 
observing projects, which have been uniformly reduced by 
ISloan et al.1 (|2003bl ) and classified based on the overall shape 
of the SED and their IR spectral fea tures (such as oxygen - 
or carbon-rich dust emission, etc.) bv lKraemer et al] (2002). 
From this database, we selected those sources with optically- 
thin, oxygen-rich dust emission (group 2.SE) sources domi- 
nated by emission from warm silicate dust (groups 3.SE and 
3. SB) and the optically-thick sources with cool silicate dust 
(groups 4.SA, 4. SB, 4.SC and 4.SE) which have a sufficient 
signal-to-noise ratio to distinguish the solid state features. 
Sources which have no significant dust emission (only SiO 
absorption at 8 pirn) , incomplete spectra and those which did 
not fulfil the criteria to be classed as an O-AGB or RSG by 
the lWoods et al] (1201 ll ) decision-tree were eliminated from 
the sample. We also require the spectra to show the amor- 
phous silicate features in order to fit their spectral energy 
distribution using the GRAMS model grid. We verified in the 

1 http : //isc . astro . Cornell . edu/~ sloan/llbrary/swsatlas/aot 1 



literature that the selected sources were classified as either 
an O-AGB or RSG. The resulting list of 110 spectra is shown 
in Table [4] In instances where a source was observed by the 
SWS on multiple occasions, we analyse each spectrum in- 
dependently. The TDT (Target Dedicated Time) number 
is used to identify the particular spectra for each source. 
Since both RSGs and O-AGB stars are red, luminous and 
exhibit similar dust emission characteristics, confusion exists 
between these two types of sources. Sources for which reli- 
able estimates for the distance are available or are known OH 
masers are further sub-divided into O-AGB and RSG stars. 
We consider distanc es from the revised Hipparcos catalogue 
(van Leeuwen|[2007l ) to be reliable if the relative error is less 
than 20 per cent. The variability types have been taken from 
the Combined General Catalogue of Var iable stars (GCVS; 
iKholopov et al.lll998l ; ISamus et al]|2004l ). 



2.3.1 Spitzer PI 094 Galactic sample 

Spitzer low-resolution spectra of 21 faint Galactic O-AGB 
stars were observed as part of Program 1094. Sp ectral data 
reduction was done in SMART (Hig don et al.ll2004l) in a man- 
ner similar to that s ummarised by lFurlan et alTl|2006l) and 
ISargent et alj l|2009h . Many of these sources are not well 
known in the literature and ancillary data including photom- 
etry, periods, variability type and distance indicators were 
difficult to find in the literature. As before we classify these 
sources using the decision tree, however, these classifications 
may be less robust than the other samples due to uncer- 
tainty in the SED shape. This means our sample may suffer 
some contamination by oxygen-rich post-AGB sources. The 
Galactic Spitzer sources are listed in Table [S] 



2.4 Globular cluster sample 

To complement these sources, Spitzer IRS observations of 
oxygen-rich sources in 14 Galactic globular clusters (GGC) 
allow the metallicity range we investigate to be extended to 
lower metallicities. The globular cluster sample consists of 39 
O-AGB stars (ob served in the Short- and Long-Low) which 
.htiflte discussed by ISloan et al.l l|2010T ) and IMcDonald et al] 
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Table 3. Details of the sources in the SMC sample. For each source the name, coordinates, Astronomical Observation Request (AOR), 
classification, bolometric magnitudes (Mbol)i period and dust mass-loss rate (Mdust) is given. 
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9 77TT 1 1 O 
Z. f / Hj — 1U 


IVLOA. OlVlVy lUO 


00 


55 


26.76 


-72 


35 


56.10 


lUOOoZoo 


ROP 


q on 
— o.UU 




1 "7RT? OO 

1. ( orj — uy 


1 1 V OOO 


00 


55 


38.30 


-73 


11 


42.00 


fiOl 7709 

uui f / yz 


ROP 
flOli 


7 91 
— / .Zl 


004 


O /lOTT 1 1 

y .4yi!j~ 1 1 


HV 11366 


00 


56 


54.80 


-72 


14 


07.00 


6018304 


O-AGB 


-6.27 


183 a 


4.23E-11 


HV 12149 


00 


58 


50.30 


-72 


18 


34.00 


6019328 


O-AGB 


-7.04 


745 a 


1.88E-09 


SMC046662 


00 


59 


35.04 


-72 


04 


06.20 


25690112 


RSG 


-8.07 


321 b 


3.10E-09 


MSX SMC 181 


01 


00 


48.18 


-72 


51 


02.10 


10665728 


RSG 


-7.28 


1062 a 


8.85E-09 


HV 11423 


01 


00 


54.80 


-71 


37 


55.00 


25691136 


RSG 


-7.78 




2.94E-09 


SMC052334 


01 


01 


54.16 


-71 


52 


18.80 


25690368 


RSG 


-7.68 




3.42E-10 


SMC055188 


01 


03 


02.38 


-72 


01 


52.90 


27603200 


RSG 


-7.54 


530° 


3.91E-09 


SMC55681 


01 


03 


12.98 


-72 


09 


26.50 


25687808 


RSG 


-7.79 




2.96E-09 


MSX SMC 234 


01 


03 


42.35 


-72 


13 


42.80 


27522048 


O-AGB 


-6.49 




4.05E-08 


HV 1963 


01 


04 


26.80 


-72 


34 


40.00 


6018048 


O-AGB 


-6.74 


330 a 


7.28E-11 


IRAS F01066-7332 


01 


08 


10.27 


-73 


15 


52.30 


17409536 


O-AGB 


-5.64 


882 c 


3.35E-08 


HV 12956 


01 


09 


02.25 


-71 


24 


10.20 


27528960 


O-AGB 


-6.27 


523 c 


5.79E-08 


MSX SMC 149 


01 


09 


38.24 


-73 


20 


02.40 


10668032 


RSG 


-7.98 




5.18E-09 


SMC083593 


01 


30 


33.92 


-73 


18 


41.90 


25690880 


O-AGB 


-6.80 




4.11E-09 



References: (al iGroenewegen et all J2009h ; (b) IPoimanskil J2002) ; (c) Soszvhski et alj J20T 



Table 4. Details of the sources in the Galactic ISO SWS sample. A full version is available in the electronic edition. For each source the 
name, coordinates, TDT (Target Dedicated Time) number, classification, distance, period and dust mass-loss rate (Mdust) ls given. 



Name 




RA 


Dec 


TDT 


Class 


Distance 


Ref. 


Period™ 


Mdust 






(J2000) 


(J2000) 






(pc) 


distance 


(d) 


(M yr" 1 ) 


S Scl 




00 15 22.18 


-32 02 43.4 


37102018 


O-AGB 


430 


16 


363 


1.35E-09 


S Scl 




00 15 22.18 


-32 02 43.4 


73500129 


O-AGB 


430 


16 


363 


1.13E-09 


WX Psc 




01 06 25.96 


12 35 53.1 


39502217 


O-AGB 


833 


1 


660 


1.92E-07 


WX Psc 




01 06 25.96 


12 35 53.1 


76101413 


O-AGB 


833 


4 


660 


1.92E-07 


OH 127.8 


fO.O 


01 33 50.60 


62 26 47.0 


44301870 


O-AGB 








2.36E-06 


OH 127.8 


fO.O 


01 33 50.60 


62 26 47.0 


78800604 


O-AGB 








6.25E-07 


SV Psc 




01 46 35.30 


19 05 04.0 


80501620 


O-AGB 


380 


11 


102 


3.08E-09 


Mira 




02 19 20.78 


-02 58 36.2 


45101201 


O-AGB 


128 


8 


332 


1.18E-08 


SU Per 




02 22 06.93 


56 36 15.1 


43306303 


RSG 


1900 


10 


533 


1.22E-08 



ta ken trom the G C JVS IKholopov et al.M199 
References: (1) Ivan Leeuwenl [1200711; (2) IChen fe Shenl [ 20081); (3) |Pe Beck et al] feoich; (4) iDecin et al.l (20071 ); 



(5) iHerman et al.1 (119861 s ): (6 ) iHumphrevsl lll978h ; (71 Jjusttanont et al.l ^OOrjl; (8) iKnapp fe Bowers! Jl98j ); (9 ) 
iLe Sidaner fe Le Bertre " 7l996ft: (101 iLevesque et alj [|2005ll; (111 lLoup et all jl993t); (121 iMassev fc Thompson! 199ll l; 
(131 lOlivier et al.l J200l1 l ; (141 lotofsson et al] l|2002l l; (151 Ivan Langevelde et al] l|l990l l: (161 IWhitelock et alj {2008 1: (171 
IXiong et alj jl994h ; (181 lYound lll995ll . 
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Table 5. Details of the sources in the Galactic Spitzer IRS sample. For each source the name, coordinates, Astronomical Observation 
Request (AOR), classification, period and dust mass-loss rate [Mdust] is given. 



Name 


RA 




Dec 


AOR 


Class" 


Period 6 






(J2000) 




(J2000) 






(d) 


(M G yr" 1 ) 


TD AC; OO^^-I-RO^I 


00 9R A1 
UU OU Zo.ftl 


UU 


A7 OQ 9^ 
ii uy.zo 


fi08981 (\ 


O-AGB 


410 
^lU.U 


K A4T7 OQ 
o . i ±'-± r_j uy 


TD AC 00 c i8Q4- l; i74'} 


01 01 A^ 


O t 


KQ A7 Sfi 

oy *± i .ou 


fiOS'3079 




1 88 
lOO.U 


K 07F 08 
O.U / IZj — uo 


TR AC 0^4.^44- ^£11 s 


C\Q A7 31 Ofi 


r >8 
■JO 


28 10 41 


fi089'}04 






1 14E— 07 


TR AC 0^14-1-9090 


0^ 1A 9fi Afi 

UO O^l ZU.'iO 


90 


99 ^ 7ti 
zz oo. t o 


aCiQ^QA 
OUOOOOft 




940 
Z^U.U 


9 QQF OQ 


TT? A 0849^ ^11 fi 


no AA r\A_ 77 
Uo H UL / I 


^1 
— 01 


97 Al 7A 


f!078Q7fl 


\ j- t \ ' 1 1 j 




4 ft^F 1 07 


TR AQ 1 Qc:Q1 ^444 


14 01 9R QO 


^4 


00 OR 
oy uu.uo 


f!081 t^fi 

uuoioou 


O-AGB 




•3 O^F 07 

O.UOHj U I 


TR AC 1 fi^9^4-.n74^ 


1 ^4 Afi A9 

It) O^i 


07 


Af) 97 07 


fi08' : i'}98 




974 
z / ^.u 


c; 71 T? OQ 
o. i I-Cj — uy 


IRAS 17030-3053 


17 06 14.07 


-30 


57 38.83 


6084864 


O-AGB 




6.58E-07 


IRAS 17276-2846 


17 30 48.30 


-28 


49 02.39 


6080768 






7.06E-07 


IRAS 17304-1933 


17 33 22.13 


-19 


35 51.75 


6084352 


O-AGB 




1.16E-07 


IRAS 17338-2140 


17 36 52.23 


-21 


42 41.08 


6085632 


O-AGB 




1.89E-07 


IRAS 17347-2319 


17 37 46.29 


-23 


20 53.48 


6080000 


O-AGB 




4.48E-07 


IRAS 17413-3531 


17 44 43.45 


-35 


32 35.28 


6081280 






1.91E-07 


IRAS 17513-3554 


17 54 42.02 


-35 


54 51.38 


6080256 


O-AGB 




3.68E-07 


IRAS 18195-2804 


18 22 40.17 


-28 


03 07.61 


6079744 


O-AGB 




4.16E-07 


IRAS 18231+0855 


18 25 33.36 


08 


56 46.62 


6082560 






1.76E-08 


IRAS 18279-2707 


18 31 03.68 


-27 


05 40.93 


6085376 


O-AGB 




2.97E-08 


IRAS 18291-2900 


18 32 19.55 


-28 


58 09.81 


6084096 






4.17E-08 


IRAS 19256+0254 


19 28 08.39 


03 


00 25.32 


6081792 


O-AGB 




1.14E-07 


IRAS 19456+1927 


19 47 49.65 


19 


35 22.67 


6081024 


O-AGB 




1.92E-07 


GX Tel 


20 11 48.22 


-55 


25 26.65 


6082048 


O-AGB 


340.8 


1.63E-08 



(a) The available data were insufficient for determining whether some sources were an O-AGB or RSG; these sources are 
not assigned a classification. 

(b) Periods are taken from the GCVS l lKholopov et alJll99St ISamus et aT1l2004h . 



(|2011bh . Both these studies compute bolometric luminosi- 
ties, periods and dust mass-loss rates for each object allow- 
ing a ready comparison to the Magellanic Cloud stars, how- 
ever, it should be noted that the mass-loss rates are obtained 
via an alternative method to that discussed in Section T3.2I 
For consistency we recalculate the mass-loss rates using the 
GRAMS models to obtain a uniform treatment of the mass- 
loss rates and directly compare the globular cluster sources 
to the other samples. Table [6] lists the sources considered in 
this work. 



3 ANALYSIS 

3.1 Bolometric magnitudes 

The bolometric magnitudes of our sources in the LMC and 
SMC samples are listed in Tables [2] and [3] These were cal- 
culated via a simple trapezoidal integration of the IRS spec- 
trum and optical/infrared photometry, to which a Wien tail 
was fitted in the optical, while a Rayleigh- Jeans tail was 
fitted to the long-wavelength data. For sources with little 
infrared excess, bolometric magnitudes we re also calculated 
using the SED-fitting code described by iMcDonald et al.l 
(2009). This code performs a \ 2 -minimisation between the 
observed SED (corrected for interstellar redd ening) and a 
grid of bt-SETTL stellar atmosphere models (jAllard et al.l 
|2011| ) which are scaled in flux to derive a bolometric lumi- 
nosity. This provided a better fit to the optical and near-IR 
photometry, than fitting a Planck function. For the most en- 
shrouded stars, fitting the SED with 'naked' stellar photo- 
sphere models leads to a underestimation of the temperature 
and luminosity due to circumstellar reddening and hence the 



integration method is preferred for very dusty sources. The 
distances and E(B-V) values adopted in our calculations are 
listed in Table [T] 

The 'classical' upper luminosity limit for AGB stars, 
base d on the core-mas s-luminosity relationship, is M;, ; = 
7.1 l|Wood et al.lll983l ); this is not an absolute limit, with 
a few AGB sour ces occasionally tr aversing this limit during 
thermal pulses (|Wood et al.lll992l ). We use this value as a 
discriminant between the O-AGB and RSGs in our sample, 
unless there is good evidence in the literature to support 
a different classification (for instance optical spectra with 
lithium absorption or OH maser emission would suggest the 
object is an O-AGB star). Fig. [JJ compares the bolometric 
magnitudes of the O-AGB and RSG sources in the globular- 
cluster, LMC and SMC samples. In the LMC we probe the 
full AGB luminosity range, however, in the SMC the ob- 
served sample is biased to the brightest AGB sources and 
supergiants. 



3.2 Determining the gas and dust mass-loss rates 

We derive individual dust mass- loss rates (Mdust) of O-AGB 
and RSG stars by fitting their infrared SED with the oxygen- 
rich subset of the 'Grid of Re d supergiant and A symptotic 
giant branch ModelS' (GRAMS; Is argent et"ail201ir ). This ap- 
proach avoids the onerous task of determining these quan- 
tities by individually modelling via radiative transfer (RT) 
the nuances of the IRS and SWS spectra. The coverage of 
the oxygen-rich GRAMS models has been compared to the 
observed colours and magnitudes of the SAGE-Spec sam- 
ple a nd was found to be in good agreement (see figs. 2-7 
from ISargent et alll201ll) . Gas mass-loss rates can be ob- 
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Table 6. Details of the sourc es in the globular cluster sample. Distances, [Fe/H] values, periods and bolometric ma gnitudes for NG C 
5139 (oj Cen) were taken from [McDonald et al. (2011b]); values for these parameters for all other clusters are from lSloan et all d201(J ) ■ 



tin g+£>t* 


**\+ q y ti *"i in o 

OLiCLl llctillL 




RA 


Dec 


AOR 


LJLo LeLllUL. 


/ 


Period. 


ivi bol 


Mdust 






(J2000) 


(J2000) 








W 




( Mr^ vr 


-l 


NGC 362 


V16 


01 


03 


15.10 


—70 


50 


32.3 


21740800 


9.25 


-1.20 


138 


-4.10 


1.84E- 


09 


NGC 362 


V2 


01 


03 


21.85 


—70 


54 


20.1 


21740800 


9.25 


-1.20 


89 


-3.50 


5.80E- 


10 


NGC 5139 


Leid 61015 


13 


25 


21.33 


—47 


36 


53.9 


27854592 


5.30 


-1.71 




-3.29 


7.92E- 


11 


NGC 5139 


Leid 42044 


13 


26 


05.35 


—47 


'V 


20.6 


27854848 


5.30 


-1.37 


0.3 


-3.01 


1.30E- 


10 


NGC 5139 


Leid 33062 


13 


26 


30.19 


—47 


24 


27.8 


27856128 


5.30 


-1.08 


110 


-3.61 


4.36E- 


09 


NGC 5139 


Leid 44262 


13 


26 


46.36 


—47 


29 


30.4 


27852800 


5.30 


-0.80 


149.4 


-3.27 


6.25E- 


10 


NGC 5139 


Leid 44262 


13 


26 


46.36 


—47 


29 


30.4 


21741056 


5.30 


-1.63 


149 


-3.88 


1.30E- 


09 


NGC 5139 


Leid 44277 


13 


26 


47.72 


—47 


29 


29.0 


27855104 


5.30 


-1.37 


124 


-3.19 


2.40E- 


10 


NGC 5139 


Leid 56087 


13 


26 


48.05 


—47 


34 


56.8 


27853824 


5.30 


-1.84 




-3.08 


1.38E- 


10 


NGC 5139 


Leid 43351 


13 


26 


55.02 


—47 


_>s 


45.9 


27854336 


5.30 


-0.98 




-2.79 


7.39E- 


11 


NGC 5139 


Leid 55114 


13 


26 


55.55 


—47 


34 


23.4 


27854080 


5.30 


-1.45 




-3.16 


1.59E- 


10 


NGC 5139 


Leid 41455 


13 


27 


15.82 


—47 


27 


54.6 


27855360 


5.30 


-1.22 


90 


-3.05 


1.64E- 


10 


NGC 5139 


Leid 35250 


13 


27 


37.73 


—47 


25 


17.5 


27855872 


5.30 


-1.06 


65 


-3.05 


4.28E- 


10 


NGC 5927 


V3 


15 


28 


00.13 


—50 


40 


24.6 


21741824 


7.76 


-0.35 


297 


-4.64 


1.10E- 


08 


NGC 5927 


VI 


15 




15.17 


—50 


38 


09.3 


21741568 


7.76 


-0.35 


202 


-4.13 


6.21E- 


09 


NGC 6356 


V5 


17 


23 


17.06 


— 17 


16 


24.5 


21743104 


15.35 


-0.50 


220 


-4.61 


1.47E- 


09 


NGC 6356 


V3 


17 


23 


33.30 


— 17 


18 


07.4 


21743104 


15.35 


-0.50 


223 


-4.06 


1.19E- 


09 


NGC 6356 


VI 


17 


23 


33.72 


— 17 


49 


14.8 


21743104 


15.35 


-0.50 


227 


-4.54 


1.33E- 


08 


NGC 6356 


V4 


17 


23 


48.00 


— 17 


18 


04.5 


21743104 


15.35 


-0.50 


211 


-4.42 


9.85E- 


10 


NGC 6352 


V5 


17 


25 


37.52 


-48 


22 


10.0 


21742848 


5.89 


-0.69 


177 


-4.27 


3.20E- 


09 


NGC 6388 


V4 


17 


35 


58.94 


-44 


43 


39.8 


21743360 


12.42 


-0.57 


253 


-4.68 


6.21E- 


09 


NGC 6388 


V3 


17 


36 


15.04 


-44 


43 


32.5 


21743360 


12.42 


-0.57 


156 


-3.83 


3.55E- 


10 


Terzan 5 


V7 




A *7 


KA QQ 
04.00 


-24 


49 


54.6 


21744128 


6.64 


-0.08 


377 


-4.98 


1.88E- 


08 


Terzan 5 


V2 


17 


47 


59.46 


-24 


47 


17.6 


21743872 


6.64 


-0.08 


217 


-4.41 


4.87E- 


09 


Terzan 5 


V5 


17 


48 


03.40 


-24 


46 


42.0 


21744128 


6.64 


-0.08 


464 


-5.09 


2.27E- 


08 


Terzan 5 


V8 


17 


48 


07.18 


-24 


46 


26.6 


21744128 


6.64 


-0.08 


261 


-5.04 


2.68E- 


09 


Terzan 5 


V6 


17 


48 


09.25 


-24 


47 


06.3 


21744128 


6.64 


-0.08 


269 


-4.72 


2.31E- 


08 


NGC 6441 


V2 


17 


50 


16.16 


-37 


02 


40.5 


21744384 


13.00 


-0.56 


145 


-3.98 


1.31E- 


09 


NGC 6441 


VI 


17 


50 


17.09 


-37 


03 


49.7 


21744384 


13.00 


-0.56 


200 


-4.26 


1.38E- 


09 


IC 1276 


V3 


18 


10 


50.79 


-07 


13 


49.1 


21742080 


5.40 


-0.69 


300 


-4.72 


1.24E- 


08 


IC 1276 


VI 


18 


10 


51.55 


-07 


10 


54.5 


21742080 


5.40 


-0.69 


222 


-4.60 


1.57E- 


09 


NGC 6637 


V4 


18 


31 


21.88 


-32 


22 


27.7 


21745664 


8.95 


-0.66 


200 


-4.44 


1.04E- 


08 


NGC 6637 


V5 


18 


31 


23.44 


-32 


20 


49.5 


21745664 


8.95 


-0.66 


198 


-4.44 


2.22E- 


09 


NGC 6712 


V7 


18 


52 


55.38 


-08 


12 


32.5 


21745920 


7.05 


-0.97 


193 


-4.34 


2.01E- 


09 


NGC 6712 


V2 


18 


53 


08.78 


-08 


41 


56.6 


21745920 


7.05 


-0.97 


109 


-3.91 


2.59E- 


09 


NGC 6760 


V3 


19 


11 


14.31 


01 01 46.6 


21746176 


8.28 


-0.48 


251 


-4.71 


6.50E- 


09 


NGC 6760 


V4 


19 


11 


15.03 


01 02 36.8 


21746432 


8.28 


-0.48 


226 


-4.34 


6.40E- 


09 


Palomar 10 


V2 


19 


17 


51.48 


18 34 12.7 


21746944 


5.92 


-0.10 


393 


-4.96 


8.85E- 


09 


NGC 6838 


VI 


19 


53 


56.10 


18 47 16.8 


21746944 


3.96 


-0.73 


179 


-4.00 


2.32E- 


09 



tained by scaling the dust mass-loss rates using an assumed 
gas-to-dust-ratio (^) based on the metallicity. Table [7] lists 
the values of ^ adopted in our calculations. We assume that 
the field stars in the Milky Way come from a single sola r 
metallicity population with ^© = 100 (Kn app et al.|[l993h . 
For the globular clusters, the gas-to-dust-ratio was deter- 
mined from the metallicity listed in Table [6] using the fol- 
lowing relation: * = 100 x l0 _ ' Fe/H l. However, it should be 
noted that the gas-to-dust-ratio metallicity dependence is 
only weakly constrained and thus the adopted gas-to-dust 
ratios are highly uncertain. 

The optical/infrared broad-band photometry for all 
sources in the Magellanic Clouds, were combined with 'syn- 
thetic' IRAS 12-fj,m fluxes obtained by convolving the spec- 
trum with the 12-/im IRAS band profilH In instances where 

2 The IRAS filte r band passes and zero points can be found in 
iNeugebauer et"aL I (QUI). 



there was more than one epoch of data for a band, the pho- 
tometry was averaged to minimise variability effects (due to 
pulsations) in the fitting. Additionally, if a source was miss- 
ing photometric data, the IRS spectrum was convolved with 
the appropriate broadband photometric filter and included 
in the SED, to better constrain the fit parameters and pre- 
vent the fit being limited by a lack of data points in a given 
region. 

The SEDs for the Milky Way sources were constructed 
using 2MASS JHK S photometry combined with IRAS 12, 
60 and 100 fj,m data, from the Point Source Catalogue 
i Beichman et alj 1 19881 ). and the Faint Source Catalogue 
I Moshirl 11989^ 7 To obtain a uniform data set we supple- 
mented the SED with convolved photometry from the ISO 
spectra in the IRAC and MIPS 24 /im bands to obtain syn- 
thetic fluxes. These were calculated using the IDL routine 
SPlTZER_SYNTHPHOT, which is available from the Spitzer Sci- 
ence Center. Synthetic photometry was also derived for de- 
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Table 7. Summary of gas-to-dust-ratios. 



Table 8. The scale factor range. 



Environment 


Gas-to-dust ratio ( W) 


MW 


100 


LMC 


200 


SMC 


500 




1585 


NGC 5139 


4266 


NGC 5927 


224 


NGC 6356 


316 


NGC 6352 


490 


NGC 6388 


372 


Terzan 5 


120 


NGC 6441 


363 


IC 1276 


490 


NGC 6637 


457 


NGC 6712 


933 


NGC 6760 


302 


Palomar 10 


126 


NGC 6838 


537 



sired photometric bands from the O-rich GRAMS models ac- 
cording to the convolution method outlined abovfl 

Variable AGB stars can show large absolute magnitudes 
amplitud es in the optical, ty pically in the range of AMy < 
11 mag l)Samus et alJ l2004T h Infrared amplitudes are less 
than this with if-band amplitudes usuall y AMk 5, 2.3 
mag and IRAC 3.6 /im AM 3 . 6 < 2.0 mag jMcQuinn et all 



120071 ; IViih et alj|200g ). Due to the large variability in opti- 
cal UBVI data we chose to exclude these points from the 
SED fitting. 

Long-period Mira variables can also show significant 
variation in amplitude due to pulsations in the mid-IR: 
multi-epoch ISO-SWS spectral studies of the variable stars 
T Cep, R Aql and R Cas covering minimum to post- 
maximum pulsation phases over a single pulsation cycle al- 
low us to characterise the variations in amplitude at longer 
wavelengths (A > 8.0 pirn). Synthetic photometry was ob- 
tained from the spectra at each epoch for the IRAS [12] 
and MIPS [24] bands, and the difference in magnitude be- 
tween photometric maximum and minimum was used as a 
measure of the amplitude. We found that AM12 < 0.4 mag 
and AM24 < 0.6. For all sources the photometric uncer- 
tainties in these bands were modified to account for this 
change in magnitude during a pulsation cycle. At near-IR 
wavelengths we also account for changes in magnitude due 
to pulsations by adjusting the photometric uncertainties to 
reflect the differences in magnitude across a pulsation cycle. 

We do not correct the SEDs for interstellar reddening. 
Due to the close proximity (d < 4 kpc) of the Galactic 
sources and the low column densities along each line-of-sight 
to the Magellanic Clouds the dereddening correction is neg- 
ligible compared to the uncertainty in the photometric mea- 
surements if variability is considered. 

A x 2-mm i m i sa tion analysis was used to fit the con- 
volved model fluxes to the observed SEDs. As the grid is 
calculated for sources at 50 kpc, we scale the convolved 
model fluxes to match those of the star before commencing 

3 In the case of 2MASS rela tive spectral response curves are pro- 
vided bv lCohen etail ||2003| 1. 



Galaxy 


Distance (kpc) 


Scale factor 


LMC 


51 ±3 


0.86 - 1.09 


SMC 


61 ±3 


0.61 - 0.74 


Milky Way 


< 4 


> 156 



the x 2 fitting. The scaling factor was calculated such that it 
minimised the difference between the source and the model, 
this allows the distance d to be a free parameter. 

The reduced y 2 , weighted by photometric uncertainties 
is given by: 



E 



(log 10 [AF^)]~log 10 [AM A (i)] 
<7(log 10 [AF A (i)]) 2 



(1) 



where \F\ is the observed broadband photometry with er- 
rors, 



a{\og m [\F x ]) = 



1 a(AF A ) 
lnlO XF X ' 



(2) 



and XM\ is the scaled model fluxes for N measurements. 
This fitting process is repeated for each model SED with a 
scaling factor and \ 2 value returned for each model. In or- 
der to investigate the onset of crystallinity determining accu- 
rate dust mass-loss rates is vital. By fitting in the log 10 [AF\] 
plane we are effectively gives more weight to the longer wave- 
length bands where the dust grains emit most of their flux 
and obtain a better fit to the shape of the SED. This method 
relies upon the source being well sampled across a range of 
photometric bands, however there are some sources where 
there is no 2MASS JHKs, IRAS [60] or MIPS [70], so the 
overall fit is more weakly constrained. 

For the LMC and SMC we assume the star lies at a 
distance d between d m i n and d max to account for variations 
in distance across each galax y due to its inclined geometry 
|van der Marel fc CionilfeoOlf ) which corresponds to a range 
of allowed luminosity scale factors for each population. For 
the globular clusters the allowed range in the scale factor 
for each source was determined, using the distances and as- 
sociated error listed in Table [5] Where possible the scale 
factor used for Galactic sources corresponds to its Hippar- 
cos distance if the relative errors are less than 20 per cent, 
otherwise we adopt a distance range of d < 4 kpc. Table [5] 
lists the distance range and scaling factor allowed for each 
galaxy. As the luminosity is a fitting parameter to the dis- 
tance, we attribute no significance to its best-fit values. 

To prevent any over-interpretation of the best fit param- 
eters, we consider all sets of model parameters that produce 



a good fit to the source, with \ 



Xbcst 



< 3, where Xbc 



is the reduced \ 2 value of the best-fitting model for each 
source , following the method employed by iRobitaille et alJ 
(2007). For each parameter we find the average value and its 
standard deviation from the range of acceptable fits. Fig. [2] 
shows an example best-fitting model to the SEDs of an O- 
AGB and a RSG in the SMC. The error bars are an in- 
dication of the expected change in flux due to variability. 
The IRS spectrum for each source is plotted to illustrate if 
a good fit has been achieved. These are not included in the 
fitting as the mineralogy is a fixed parameter in the GRAMS 
model grid. 
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Table 9. Wavelength intervals for extracting crystalline silicate 
emission features cf . Fig. [4] 



Wavelength [micron] 

Figure 2. Best-fitting models to the broadband photometry for 
a RSG (MSX SMC 234; dashed red) and an O-AGB star (MSX 
SMC 149; dashed blue) in the SMC. The solid line is the observed 
spectra for each star. Similarly the square and diamond points 
are the respective observed photometry for the O-AGB and RSG 
sources. The photometric errors have been estimated based on 
the variability. 



The dust mass-loss rates determined through the spec- 
tral fitting of the SEDs with the GRAMS model grid are listed 
in Tables [2HS1 In principal, the fitting routine also provides 
uncertainties for the individual dust mass-loss rates, of or- 
der 10 to 55 per cent, however, in absolute terms, the errors 
are much greater than this due to the inherent uncertainties 
and assumptions in the model. 

In order to assess the accuracy of the absolute values 
derived in the GRAMS fitting, we comp are our derived dust 
mass- loss r ates to those dete r mined by iGroenewegen et all 
(|2009l ) and iMcDonald et al.l (|2011bl ) for common sources 
between the data sets. The GRAMS fits are systematically 
greater by a factor of 3-8. This is predominantly due to 
a difference in the choice of optical constants and grain 
sh ape. The gram s grid uses the as tronomical silicates 
of Osscnkopf et al. ( 1992jV while both iGroenewegen et al.l 
(|2009l ) and lMcDonald et al l (|2011bl ) use a combination of sil- 
icates, metallic iron and aluminium oxide. Different physical 
parameters used in modelling the dust envelope and central 
star also contribute to this discrepancy. For more details on 
how the GRAMS mod els compare to other modeling efforts 
we refer the reader to lSargent et al.l l|201ll ). 

Our analysis regarding the onset of crystallinity is very 
sensitive to the uncertainties in the dust and gas mass-loss 
rates, thus it is important to account for the possible sources 
of error and caveats used in the modeling to prevent misin- 
terpretation. In the models we have assumed that all out- 
flows have the same constant velocity of 10 kms -1 , whereas 



Feature (pm) 


Continuum interval (/im) 


23.6 


22.3 - 24.6 


28.1 


26.0 - 31.1 


33.6 


31.5 - 35.0 



in reality they may have expansion velocities between ~5 
and ~25 kms - , even within the same galaxy. However, 
by assuming a constant velocity the dust and gas mass-loss 
rates obtained in the fitting can be used as a proxy for the 
gas and dust column densities. As we have been consistent 
in our treatment of the samples and in the determination 
of the dust mass-loss rates, the relative uncertainty between 
the sources is small compared to the absolute uncertainties, 
allowing comparisons to be drawn across the sample. The 
derived gas mass-loss rates have an additional uncertainty 
because of the adopted scaling relation between dust and 
gas in circumstellar environments. For oxygen-rich evolved 
stars, this relation depend strongly on metallicity, however, 
it is not well constrained. As such, the absolute values for 
the gas mass-loss rate should be treated with a degree of 
caution, nevertheless, we expect the assumed gas-to-dust ra- 
tios to provide a reasonable first order approximation of the 
relative gas mass-loss rates. 



3.3 Measuring the crystalline silicate features 

Narrow spectral features due to crystalline silicates tend to 
cluster in well-defi ned complexes, found near 10, 18, 23, 28, 
33, 40 and 60 (im dMolster et al.ll2002bf ). We do not extract 
information for the features near 10 and 18 fim as they are 
impossible to cleanly separate from the broader, underlying 
amorphous silicate features. 

In order to determine the properties of crystalline sil- 
icate emission features across a range of metallicities, we 
have defined a continuum representing a featureless ther- 
mal dust component, which we used as a baseline to obtain 
an enhanced view of the spectral features beyond 20 fim. 
To produce a smooth continuum, a 5-pixel boxcar smooth- 
ing algorithm was applied to the IRS spectra to reduce the 
noise fluctuations between each data point (while preserving 
the resolution). However, it is possible that this may result 
in any low-contrast, broad dust features being treated as 
continuum. The continuum was determined by applying a 
cubic spline fit in F„-space for A > 19.5 /im. Spline points 
were fitted to the spectra at 19.8, 21.3, 25.5, 31.3, 35.6 /im 
where the continuum is well defined and no strong spectral 
features are present. Examples of the spline fitting proce- 
dure can be seen in Fig. [3] In instances where a feature is 
present in a spectrum at a standard spline point (for ex- 
ample SSID130), the spline point was moved to the closest 
alternative wavelength where no feature could be visually 
identified. A similar procedure was applied to the ISO SWS 
spectra after regridding to the Spitzer IRS resolution. 

We prefer the above method compared to the alternate 
approach used by some authors, in which a local continuum 
is defined for each feature via linear interpolation between 
flux densities in a given wavelength range on either side 
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1.08 




20 25 30 35 

Wavelength [»m] 



Figure 3. Example of the spline fitting for an LMC O-AGB 
(top), LMC RSG (2nd panel), SMC O-AGB (3rd panel), SMC 
RSG (4th panel). The spline (dashed line) is fitted to the spectra 
at the marked wavelengths. The top three spectra show crystalline 
emission features, while the spectrum in the 4th panel does not. 
The bottom panel shows SSTISAGEMC J053128.44-701027.1 
(SSID 130), where alternative spline points were selected in-order 
to produce a smooth continuum. Fig. [4] presents the resulting 
continuum-divided spectra. 

of the feature. As slight variations in the position of the 
feature and noise in the spectra all have strong effects on the 
determination of the continuum because the fitting region 
is small. Also the dust continuum is inherently non-linear, 
with changes in curvature, producing a non-systematic offset 
between sources. 

We determine equivalent widths of the 23-, 28- and 33- 
/im features in each spectrum to provide a quantitative mea- 
sure of their strengths. These are calculated with respect to 
the underlying continuum as determined by the aforemen- 
tioned spline fitting process. The equivalent widths (W cq ) 
can then be determined by: 

W eq = Wl-|^)dA, (3) 

where F ^ s and -Fcont are the observed flux and the flux of 
the spline-fitted continuum. The wavelengths used to define 
the blue and red sides of each feature are given in Table 
[9] An example of the feature extraction process is given in 
Fig.H 
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Figure 4. Continuum-divided spectra from the sources shown in 
Fig. [3] The dotted vertical lines mark the regions within which 
the complex strength was measured. 

If a feature is particularly weak compared to the contin- 
uum, its strength may be influenced by how the continuum 
is fit. To quantify this effect, we also determined equivalent 
widths using linear interpolation to approximate the con- 
tinuum under the features. For spectra with high contrast 
features or with a good S/N ratio the difference in strength 
is of the order of < 15 per cent, however, there is a large 
disparity (up to ~65 per cent) in strengths for low contrast 
features. 

Both the spline and linear interpolation methods will 
underestimate the total contribution from the crystalline sil- 
icate components to the total flux. These methods fit the 
continuum under each crystalline silicate emission feature, 
however, the crystalline silicates grains will also contribute 
some flux to the underlying continuum which is composed 
of contributions from all dust species of varying mineralogy, 
lattice structure, temperature and size and shape distribu- 
tions. 

4 RESULTS 

4.1 Mass-loss rates 

We have separated the O-AGB and RSGs stars in each 
galaxy in order to examine the influence of the mass-loss 
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Table 10. Summary of mean mass-loss rates. 
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(a) M^ust distribution (b) M distribution 

Figure 5. Overview of the dust and gas mass-loss rate distributions for each population. Column 1: O-AGBs with solely amorphous 
silicates, column 2: O-AGBs with crystalline silicates, column 3: RSGs with only amorphous silicates, column 4: RSGs with crystalline 
silicates. The box encloses the inter-quartile range, the whiskers extend out to 1.5 times the lower or upper quartile or to the maximum 
or minimum value of the data (if smaller). Outliers are indicated by circles. 
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Figure 6. Bolometric luminosity versus dust mass-loss rate for 
the LMC (red colours), SMC (blue) and globular cluster (green) 
sources. O-AGB stars are plotted as triangles and RSG as circles. 
Objects with crystalline silicates are plotted with closed symbols. 
Milky Way sources are not included because of distance uncer- 
tainties. 



process on the dust production in these two classes of object 
separately. Furthermore, we distinguish between sources in 
which crystalline silicates are visually apparent and those 
whose silicates appear to be solely amorphous. For each 
population we determine the mean log Mdust and logM. In 
general, the dust mass-loss rate for the sources with only 
amorphous silicates increases as the metallicity of the sam- 
ple increases. Thus, the SMC and globular cluster sources 
have low dust mass-loss rates while the Galactic field sources 
have the largest dust mass-loss rates. The individual objects 
within a galaxy encompass a range of mass-loss rates as in- 
dicated in Fig. [5] An overview of (log Mdust) and (log M) 
for each population is presented in Table [TO] 

In each instance the mean dust and gas mass-loss rate is 
higher for the sources which exhibit crystalline silicate fea- 
tures compared to those that do not. Also, while a low mass- 
loss rate is no guarantee that the silicates will be entirely 
amorphous, it is evident that high mass-loss rate sources 
have a greater disposition to exhibit crystalline silicate fea- 
tures. The mean dust mass-loss rates for the sources with 
crystalline silicates do not show the same metallicity depen- 
dence as the sources with only amorphous silicates. As these 
objects have different gas densities but similar dust column 
densities this may indicate that crystallinity depends on the 
dust density in the wind. The mass of the central star may 
also be significant: a greater fraction of O-AGB stars (56 per 
cent) display crystalline bands compared to the RSGs (24 
per cent) and, as the metallicity decreases, the proportion of 
RSGs with crystalline features also drops. This may be due 
to differences in the silicate dust composition of RSG and 
O-AGB stars, with the former exhibiting amorp hous Ca-Al- 
rich s ilicates before formin g Mg-rich silicates |Speck et al.l 
l200d ; IVerhoelst et aljlioosh . 

Fig. [5] shows the bolometric magnitude as a function of 
dust mass-loss-rate. We detect crystalline silicates across the 
full AGB and RSG luminosity range, from the early AGB 



to the most luminous RSGs. This indicates that we are sen- 
sitive to crystalline features down to a low luminosity limit, 
not just in the brightest sources where there is better con- 
trast between the feature and the noise. For the Magellanic 
Clouds our sample is biased towards the brightest AGBs 
and RSGs or dusty sources with high mass-loss rates. Here 
the threshold density for the onset of crystallinity is greater 
than that of the Galactic field stars and globular cluster 
sources. As such, the detection of crystalline features in the 
Magellanic Clouds may be limited by the signal-to-noise. 



4.2 Mass-loss rates and crystalline feature 
strength 

Fig. [7] shows the strengths of the 23-, 28- and 33-/im fea- 
tures against the dust mass-loss rate (upper panel) and total 
mass- loss rate (lower panel). In both instances, objects with 
crystalline silicate features are represented by filled symbols, 
while open ones represent sources where crystalline silicates 
are not visually apparent. It is interesting to compare and 
contrast the behaviour of the crystalline features at the dif- 
ferent wavelengths, in particular the 23- and 33- /im features, 
which have the same principal carrier, forsterite. At 33 /im, 
features due to the the crystalline silicate grains only become 
apparent at the highest mass-loss rates, where the major- 
ity of the sources display some crystallinity. This is usually 
interpreted as evidence that there is a threshold value in 
the mass-loss rat e at which crystalline silicates can form 
|Cami et alj|l998t ISvlvester et al.lll999h . We observe, how- 
ever, that the forsterite feature at 23 /im paints a different 
picture. Here, the transition between sources which exhibit 
a feature and those which do not is blurred, indeed at the 
highest mass- loss rates (M > 10 -4 MQyr -1 ) the detection 
rate of the 23-/tm feature drops: this may be a consequence 
of the forsterite grain temperature or it may indicate that 
the 23-/im band is going into absorption due to high column 
density. An example of the changing relative strengths of 
the two forsterite features as the optical depth of the shell 
increases is plotted for three MW stars in Fig. [8] Like the 
23-/tm forsterite feature, the 28-/xm enstatite feature spans 
a wider mass-loss rate range than the 33- ^m feature. How- 
ever, in all three cases no evidence for crystallinity is seen 
in the lowest mass-loss rate sources. 

The sources where no feature is visually apparent 
are well described by a Gaussian distribution around a 
strength/continuum ratio of zero. This implies that this 
spread is due to stochastic noise in the spectra and that we 
can trust systematics in the crystalline silicate detections. 
For the 23- and 33-//m bands there is a slight positive skew, 
this is probably due to real features which we are unable to 
visually identify in the data as a consequence of substantial 
noise in the Spitzer spectra at these wavelengths. The 33- 
/im band can be susceptible to contamination from diffuse 
[S III] and [Si II] emission lines, which may be improperly 
subtracted due to their non-uniform variation across the slit. 
We have taken care to minimise this contamination by ex- 
cluding this region from our analysis if we detect the pres- 
ence of sky lines. Consequently, the detection rate of the 
complexes represents a lower limit. 
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Figure 7. Top: The strength of the 23-, 28- and 33-^tm features versus Ma us t, the dust mass- loss rate. O-AGB stars are represented 
by triangles, RSG as circles and sources which cannot be distinguished as either O-AGB or RSG are indicated by squares. The LMC 
sources are plotted in red; SMC: blue; MW: orange; and the globular clusters: purple. Filled symbols indicate that a crystalline feature 
can be visually identified. For sources where no feature is detected, a Gaussian is fitted to the distribution of feature strength/continuum. 
Bottom: The strength of the 23-, 28- and 33-fj.m features versus M, the total mass-loss rate. As with the dust mass-loss rates there is no 
apparent threshold value between sources with and without 23- and 28-fim silicate features. 



4.3 Dust temperature 

The strengths of individual bands are affected by the abun- 
dance of each species, and by grain shape and temperature. 
By eliminating the former two parameters, we can derive the 
temperature of the emitting grains. We can eliminate the 
need to parametrise the abundance by taking the strength 



ratio of features which have the same carrier: in the case of 
forsterite, we can use the ratio of the 23- and 33-/im features. 

Eliminating the effects of grain shape are more difficult, 
as these are further confounded by the choice of laboratory 
optical constants used to create the feature shapes. Emission 
models for the forsterite gra i ns were computed from labora- 
tory data (jjager et al.lll998l ; iKoike et al.fll99Sl ). Both these 
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Figure 8. Three Spitzer spectra ordered in terms of the 10 
optical depth. As the dust shell becomes optically thick the 23-/im 
feature strength becomes suppressed. The smooth red lines show 
the adopted continua (spline). The shaded regions indicate the 
wavelengths where the 23- and 33-^tm features were measured. 
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Figure 9. The ratio of relative peak strengths above con- 
tinuum for the 23- and 33- fim forsterite bands as a function 
of the forsterite dust grain temperature. The different curves 
represent the diffe rent forsterite laboratory data: s olid purple: 
l|Koike et al] H998T ): dashed red: <|jager et al.lll998h continuous 
distrib ution of ell ipsoids shape distribution (CDE); and dash dot 
blue: l|Jager et allll99Sh homogeneous spherical particle shape 
distribution (Mic). 



datasets give a good match to the observed feature profiles 
and do not differentiate different crystallographic axes, re- 
ducing the number of parameters which we must fit. From 
these, we use three sets of a bsorption efficienci es (Qf rs(^))- 



The first comes directly from lKoike et al 



dl998l). The secon d 



Jager et all |l998), 



is derived from the optical constants of 
using homogeneous spherical particles with radii of 0.01 fim 
via Mie theory. The third is similarly calculated using a con- 
tinuous distribution of ellipsoid grains with equivalent vol- 
umes. 

Assuming that the long wavelength (A > 21 pun) 
forsterite features are formed in an optically-thin region of 
the outflow then the crystalline forsterite spectrum can be 
approximated by: 



F(u) 



B(T ioTB ,is) x Q(y) 

fors 



(4) 



where B(T{ OIS ,iy) is a blackbody of temperature Tf ors , 
QMfors is the forsterite absorption efficiency, and F con t 
represents the sum of the featureless thermal dust compo- 
nent and stellar continuum. This component is effectively 
removed from the spectra by subtracting the spline-fitted 
continuum we use in Section 13.31 We therefore also subtract 
an identically-fitted spline from the modelled forsterite spec- 
trum to ensure consistency. 

These continuum-subtracted forsterite absorption effi- 
ciencies were then multiplied by blackbody functions corre- 
sponding to temperatures between 30 and 1000 K at 5 K in- 
tervals. The strengths of each feature were then determined 
by summing the remaining flux over the wavelength inter- 
vals given in Section ^, 31 and a ratio of the 23- to 33-/im fea- 
tures was determined. Fig. [9] shows how this ratio changes 
with temperature, allowing direct comparison to observed 
band strengths. An average of the resulting forsterite tem- 
peratures are given in Table 1111 In sources where there is 
no 23-ptm feature detection, T/ ors may be incorrect due to 
self-absorption. This approach relies on the ratio of features' 
strengths, it is not greatly affected by the change in feature 
position due to temperature, and can therefore be treated as 
a robust estimate of temperature once one takes into account 
the above caveats regarding grain shape and self-absorption. 

In the 23-/xm wavelength regime, the suppression of the 
forsterite feature for the majority of the sources where crys- 
tallinity is present is probably a consequence of the poor 
contrast of the cool grains above the amorphous dust con- 
tinuum and, as such, they are indiscernible with Spitzer and 
ISO. For the sources with extremely high mass-loss rates 
we believe the absence of the 23-pim forsterite feature is 
due to self-absorption. Fig. QJj] shows the spectra of four 
Galactic O-AGB stars undergoing intense mass- loss rates. 
In these sources, amorphous silicate absorption features are 
seen at 10 and 20 ^m; OH 32.8-0.3, OH 26.5+0.6 and 
IRAS 17276—2846 also clearly exhibit crystalline silicate ab- 
sorption features up to 33 fj,m where the feature appears in 
emission. The extreme OH/IR star IRAS 16342-3814 is also 
know n to have some crys talline silicate features in absorp- 
tion l|Diikstra et al.| [2003). This object is undergoing an in- 
tense mass-loss of ~ 10 -3 Mq yr _1 , which is approximately 
an order of magnitude greater than the highest mass-loss 
rates determined in this work. This be haviour is similar to 
that modeled by Ide Vries et al.l l|2010l ) , who found that at 
mass-loss rates greater than 10 -5 Mq yr _1 even the 33-pim 
band becomes self-absorbed in some lines of sight. These 
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Figure 10. Spectra showing crystalline silicate absorption and 
emission features. To enhance the visibility of the crystalline fea- 
tures in the spectra the adopted continua are shown (smooth red 
lines). The position of the 23-^tm feature is indicated by a dashed 
line. 

results show that there is a limit to the strength of the crys- 
talline silicates in AGB stars. 

4.4 Feature shapes & positions 

To compare the crystalline features across the Milky Way, 
Magellanic Clouds and the Galactic globular clusters we de- 
rive a mean continuum-subtracted spectra for the complexes 
at 23, 28 and 33 fim. The emission from the crystalline sil- 
icate features can be isolated from the observed spectrum 
by subtracting a continuum according to Equation [4] Emis- 
sion features were extracted from the resulting continuum- 
subtracted spectra over the wavelength ranges specified in 
Table [9] and averaged using a weighting factor (calculated 
from the point-to-point RMS) to account for the varying 
signal-to-noise in each region. The mean spectral complexes 
for the MW, LMC and SMC sources are shown in Fig. [11] In 
the globular cluster sample only the 28-^m spectral complex 
is present. 

The precise shape and position of the continuum- 
subtracted features will depend on the spline-fitted contin- 
uum. In order to evaluate how sensitive a feature is to the 
defined continuum, we also computed average profiles using 
a linear continuum. The central wavelength of the 23- /im fea- 
ture (measured by finding the wavelength which bisects the 
integrated flux of the feature) has a mean variation of 0.08 



Table 11. Derived forsterite temperatures for stars with a 33- 
fim feature. In sources where there is no 23-fim feature detection, 
Tf ors may be incorrect due to self-absorption. 



Star 23 


fim Ft. 


T fors (K) 


NMT, Pvcr 


Y 


88 zb 10 




Y 


121 zb 23 


WX Psc 


Y 


89 zb 10 


WX Psc 


Y 


89 zb 10 


TRP 4-^01 37 


Y 


170 zb 47 




Y 


102 zb 13 


OH 1 04 Ql 4-2 41 


Y 


67 ± 6 


TR AC 034344-^81 8 


Y 


luO zrz ouu 


TR AC, 04^3 fiR9^ 


Y 


73 4- 7 
i o zrz f 


TRAS 04^4^—7000 


Y 


88 zb 10 


TRA , 0^409— fiO^ 


Y 


98 zb 12 


TRA^ 0^908— ^0^7 


Y 


108 zb 15 


TRAS 0^380—^1099 


Y 


4QC 4_ Qqn 
iyo zrz oyu 


TRAS 13^811—^444 


Y 


1 QQ 4_ ftQ 

iyo zn uo 


TRAS 17304—1033 


Y 


QOI 4_ 1Q0 

o^-L zrz iyu 


TRAS 17338—9140 


Y 


1 CQ 4_ Af, 

iuy zrz iu 


TRAS 17413—3^31 


Y 


1 rro _|_ on 
too zrz oy 


TRAS 18970—9707 


Y 


r:r)7 4_ ooc 
ou i zrz ooo 


TR AC 18901 9000 


Y 


94^ 4- 108 
z^o zrz iuo 


TRAS 109^^4-09^4 


Y 


r\Af\ 4- 3R4 
o £ ±u zrz oui 


TR AC 1QA^fi4_1Q97 


Y 


o«7 4- 1 30 
zu i zrz iou 


M^V ^.A/TP 1 34 


Y 


117 4_ 10 

iii zrz iy 


MSX T,MP 807 


t\ 


< 150 


S Scl 


\ 


< 64 


HD 9fi0^0Q 




< g»g 


ApnT, 91 00 
ill vjtj ziyy 




C 70 


ATTPT 9403 




< ^3 


AT?nT, 9403 


\ 


<; Qg 


A PPT 9374 




< 8Q 

5^ oy 


ATTPT ^370 


fs; 


< A7 


ATTPT ^370 


fs; 


<; 


A PPT ^3^ 


f^- 


<T ^0 


OH 91 ^4-0 ^ 


f\ 


C 52 


OH 9fi 9—0 


]\" 


<^ 43 


OH 9fi ^4-0 R 


\ 


<; Q3 


OH 39 8—0 3 
v./ 1 i oz.o — yj.o 


t\ 


< 45 


OH 1 97 8 4-0 


N 


^ 84 


OH 127.8 +0.0 


N 


^ 56 


IRAS 05128-6455 


N 


^ 70 


IRAS 04407-7000 


N 


^ 122 


IRAS 08425-5116 


K 


^ 56 


IRAS 05329-6708 


N 


^ 206 


IRAS 05558-7000 


N 


^ 153 


IRAS 17004-4119 


N 


^ 56 


IRAS 17010-3840 


N 


C 51 


IRAS 17030-3053 


N 


C 53 


IRAS 17276-2846 


K 


^ 54 


IRAS 17347-2319 


N 


^ 47 


IRAS 17513-3554 


N 


^ 84 


IRAS 18195-2804 


N 


C 60 


IRAS 18231+0855 


N 


^ 513 


MSX SMC 024 


N 


C 84 



fixa when a linear continuum was fitted, a mean change in 
position of 0.19 /im was recorded for the 28-/im feature and 
the 33-/im feature had a 0.17 /im change. In extreme cases 
where the features are weak compared to the continuum, 
or the spectra have a low signal-to-noise the central wave- 
length may move by 0.43 /zm. For most cases the change in 
wavelength is relatively minor. The average feature shape is 
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Figure 11. Overview of the normalised mean continuum-subtracted 23-, 28- and 33-/im complexes. Top: the LMC, SMC, globular 
cluster and MW Spitzer complexes are plotted in blue, red, green and brown respectively, while the re-gridded MW spectra from the 
SWS are plotted in orange. Middle row: the combined mean complex spectrum of O-AGB (blue) and RSG (Red) sources. The tick marks 
indicate the wavelengths where crystalline features are found. The crystalline species are indicated by a F (forsterite) , E (enstatite), D 
(diopside), whil e a U indicates an unidentified species. Bottom: Laboratory spectra of forsterite (black) and enstatite (grey) are plotted 
for comparison (Koike et al. 1 9981 ), 



remarkably insensitive to the method of continuum fitting, 
due in large part to the weighting factor and normalisation. 



4-4-1 The 23-fj.m complex 

The 23-^m complex is clearly seen in several of the 
stars across a wide range of dust mass-loss rates. How- 
ever, it becomes less prominent and even absent in the 
highly-enshrouded sources. This band is entirely absent in 
the globular cluster sources. The co mplex is thought to 
be dominated by forsterite emission (|Molster et all [2002b; 
ISloan et al.ll2008r ). with its peak fitting the dust spectra bet- 
ter than Fe-rich members of the olivine series. The average 
central wavelength of the complex, for each population, is 
around 23.35 /im, however, the variations between samples 



and individual sources is larger than that of the 33-/xm fea- 
ture. 

The complex is comprised of several blended features at 
22.3, 22.8, 23.0 and 23.7 /im whose relative strengths provide 
the main difference between the samples. While the 22.3, 
23.0 and 23.7 i xm feat ures have previously been reported by 
iMolster etafl (|2002bh . this is the first-time the 22.8 (im com- 
ponent is seen. The complex is dominated by the 23.7-/xm 
forsterite band, with the 23.0-/im enstatite feature also con- 
tributing to the substructure. The 22.3-/mi and the newly 
discovered 22.8 /im bands have yet to be attributed to a dust 
species. 

The strength of the 23.0-/xm feature compared to the 
23.7-/im feature provides the main difference between the 
Magellanic Cloud mean profiles and the MW profiles. The 
mean profile of the SWS Galactic sources exhibits the largest 
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deviation from the other mean spectra: this mean spectrum 
shows prominent 22.3- and 22. 8- rim features which are more 
subdued in the other profiles. Furthermore, the mean SWS 
Galactic spectrum is shifted towards the blue. The mean 
profile of the RSGs has a similar shift to shorter wavelengths 
with a central wavelength of ~23.0 fim (compared to the O- 
AGB sources at 23.4 /-tm), as seen in the middle left panel 
of Fig. [Til 

4-4-2 The 28-fim complex 

The 28-/im complex shows clear variations in the spectra 
with metallicity. Although the samples show a similar peak 
position around 27.8 /jm, the spectral shape is different for 
the Magellanic Cloud, globular cluster and Galactic sources. 
This complex is comprised of several blended features of 
which the 27.5-iim forsterite feature, combined with the 
28.2- and 29.4-/xm enstatite features, are the main compo- 
nents |Molster et al.ll2002bl) . In the Magellanic Cloud pro- 
files, the broader feature may indicate an enhanced enstatite 
contribution compared to forsterite. There may also be a sig- 
nificant contribution from the unidentified 29.6 /im feature, 
which is especially pronounced in the MW Spitzer mean 
complex. Unlike the MW and the Magellanic Cloud pro- 
files, the mean spectrum for the globular clusters does not 
show a forsterite component; here, the complex is dominated 
by enstatite. There is no discernible difference between the 
O-AGB stars and the RSGs. 

Determining an average profile for the MW SWS spec- 
tra was hampered by the troublesome boundary between 
bands 3D, 3E and 4 in this region, where known light 
leaks and spectral response cali bration prob l ems re sult in 
large uncertainties in the flux (|Sloan et al.l l2003 aO . Sev- 
eral SWS spectra show a broad plateau extending from 
~ 27.5—32.2 ixm. We could find no evidence for a similar- 
feature shape in any of the Spitzer spectra and we conclude 
that this plateau is probably an instrumental artifact. We 
have taken care to exclude any sources where this artifact 
may be present in the ISO SWS data before determining a 
mean spectrum for this sample. 

There is a large discrepancy between thes e mean profiles 
and the 28 /xm outflow complex determined bv lMolster et al.l 
(2002b). That complex shows a sharp rise across the band, 
furthermore the red-edge of the complex has a strong rel- 
atively flat component with few peaks, which is not seen 
in any of the Spitzer spectra. We believe that our 28-^tm 
profiles provide a more accurate/reasonable reflection of the 
crystalline features present, as these are not affected by the 
problematic SWS 3D and 3E bands. 



4-4-3 The 33-[im complex 

The 33-pm complex is dominated by the strong forsterite 
feature at 33.6 \ivn, and is particularly prominent in sources 
with high mass-loss rates. In both the LMC and SMC sam- 
ples, strong noise around 33 (im complicates the identifi- 
cation and extraction of features, thus the detection rate 
may be lower than expected compared to the MW sam- 
ple. The 33-/im complex is not seen in the globular cluster 
stars. The complex shows remarkably little variance with 
metallicity, both in terms of the shape, mean wavelength 




Wavelength [um] 

Figure 12. Spitzer spectra of three sources which display very 
strong crystalline features. We postulate that these sources have 
a stable disc-like geometry. 

and FWHM. The mean wavelength of the extracted feature 
is 33.54 (im, varying by 0.13 /im across the three galax- 
ies. The position of the 33.6-/im feature is consistent with 
cool Mg-dominated crystalline olivine grains. If the grains 
contained iron incl usions the peak of the feature would b e 
shifted to the red (|Koike et al.l 120031 ; iPitman et al] 120101 ). 
Warmer grains would also result in a shi ft in band position 
to longer wavelengths (|Koike et al.ll2006h . The 33-/xm com- 
plex also exhibits some less prominent sub-structure. There 
is a weak feature on the left shoulder of the complex at 
32.2 ^m, tent a tively identified as diopside (MgCaSi2 0e) by 
iMolster et all (|2002bl ). At the red edge of the band we note 
the presence of an enstatite feature at 35.3 /xm. 

4-4-4 Probable disc sources 

Fig. shows three sources that exhibit unusually strong 
crystalline silicate features in their spectra. Strong crys- 
talline silicate features such as these are often found in 
long-lived circumstellar discs where dust pro cessing creates 
very high abundances of crys talline silicates l|Molster et alj 
ll999l ; lGielen et al.ll2008l , l201ll ). We propose that HD 269599 
(SSID 4464), MSX SMC 134 and HV 12956 have a sta- 
ble disc-like geometry, which invalidates the assumption of 
a spherically-symmetric outflow for the SED fitting. HD 
26959 9 is a member of th e young (16 ± 2 Myr) cluster NGC 
1994 (Ku mar et al.l [20081 ) which contains B and B[e] stars. 
Because of its large infrared excess it is probable that HD 
269599 is also a relatively massive, relatively young star with 
a disc. Optical spectroscopy confirms that H D 269599 is a 
B[e] supergiant (|Lamers et al.|[l99 8l; lKastner et al.ll20ld ). 
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Figure 13. Correlation between the strength of the 23-fim fea- 
ture and the 28-fj.m feature. The symbols are as defined in Fig. [7] 
except the filled symbols represent the presence of both a 23- and 
a 28-fim feature. The lines show the least square fit to the LMC 
(red) and MW (orange) features. 




0.01 0.10 

33 feature/continuum 

Figure 14. Correlation between the strength of the 33-fim fea- 
ture and the 28-fim feature, with best-fit line for each population. 
The symbols are as defined in Fig. [7] except the filled symbols 
represent the presence of both a 33- and a 28-fim feature. 



MSX SMC 134 is RAW 631, an optically identified car- 
bon star. Optical and near-infrared spect roscopy confirm 
its ca rbon-rich nature (Sloan et al. in prep, Ivan Loon et al.l 
2008). The combination of a carbon-rich photosphere and 
oxygen-rich dust is the hallmark of a si licate carbon star 
(first discovered by iLittle-Mareninl 1 19861 ) . The leading ex- 
planation for these sources is that the silicate dust was 
trapped in a disc in a binary system before the mass- 
losing star evolved in to a carbon star l|Llovd Evans! Il99ll ; 
iBarnbaum et al.lll99ll ). 

We could find nothing in the literature that signifies 
the presence of a disc around HV 12956, but postulate the 
presence of a disk is required to show these high-contrast 
crystalline silicate features. 

4.5 The olivine-to-pyroxene ratio 

The variation in spectral shape and the differences between 
individual bands present in the 28-pim complexes for the 
Milky Way, Magellanic Clouds and globular clusters, may 
indicate a change in the crystalline silicate dust mineralogy 
with metallicity. As the metallicity increases, the enstatite 
substructure at 28.2 and 29.4 /im becomes less pronounced, 
whereas the forsterite substructure at 27.5 /im becomes more 
distinct. This dichotomy between the Magellanic Clouds and 
the Galaxy can be clearly seen in Figs.[l3]and[l4]which show 
the relations between the 23- and 33-/im forsterite features 
and the 28-/im enstatite feature strength. The separation 
between the populations, as indicated by the difference be- 
tween the lines of best fit to the crystalline silicate features 
in each galaxy, suggests a change in crystalline dust compo- 
sition with metallicity. A line of best fit is not plotted for 
the SMC sources due to insufficient data points. Note that 
the globular cluster sources do not feature in these plots 
as they only exhibit the 28-pm enstatite feature, also indi- 
cating a change in dust composition with metallicity. This 




0.01 0.10 

33 feature/continuum 



Figure 15. Correlation between the strength of the 23-fim fea- 
ture and the 33-fim feature. The symbols are as defined in Fig. [7] 
except the filled symbols represent the presence of both a 23- and 
a 33-fj.m feature. 



trend is almost certainly a metallicity effect, and not due to 
changes in mass-loss rate or luminosity as we are sensitive to 
crystalline silicate features in AGB and RSG sources across 
several dex in both luminosity and mass- loss rate. 

Figure [T5l shows the relation between the 23- and 33- ^m 
features which have the same primary carrier (forsterite). 
There is no clear separation between the different metallicity 
populations. However, several of the Milky Way sources are 
not very well represented by the best fit line. These sources 
are experiencing intense mass-loss rates consequently the 23- 
/j,m feature might experiencing self-absorption. 
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Figure 16. Comparison between pulsation period and the crystalline silicate feature strength. Symbols are the same as Fig. [7] 



4.6 Pulsation period and crystalline feature 
strength 

In general, as an AGB star evolves, the stellar radius 
increases, it becomes brighter, and its pulsation period 
increases. A similar period-luminosity relation exists for 
RSGs, although they often have smaller amplitudes and less 
regular periodicities. Fig. \W\ compares the crystalline band 
strength with pulsation period. Care needs to be taken when 
comparing trends with pulsations across a metallicity range, 
since a star's effective temperature is metallici ty-dependen t 
and hence so are radius and pulsation period |Woodlll990l) . 
We therefore consider each population separately. We find 
that the fraction of sources showing crystalline silicate fea- 
tures is greater for sources with pulsation period longer than 
700 days than for sources with short pulsation periods. Fur- 
thermore, in these long period pulsators the strength of the 
crystalline features appears to be increasing with period. 
This effect is strongest for the 28 /im feature in Fig. 1161 
These stars are typically more evolved with higher dust 
mass-loss rates than the sources with periods less than 700 
days. 



5 DISCUSSION 
5.1 Mass-loss rates 

Determining total mass-loss rates requires an assumption 
about the scaling relation between the dust and gas in the 
envelope of AGB stars; this correction factor is not well 
known, even fo r Milky Way sources. It was proposed by 
Ivan Loonl (2000) that the gas-to-dust ratio for O-rich sources 
is approximately linearly dependant on the initial metal- 
licity. Although this relation has an inherent scatter, it is 
commonly applied when deriving gas mass-loss rates. Our 
estimates for the total mass-loss rates are similarly affected 
by the large uncertainties in the correction factor, which 
may blur the effects of the gas density on the onset of crys- 
tallinity in the envelopes of evolved stars. Furthermore, the 



precision of the individual gas mass-loss rates will be af- 
fected by the assumption that objects in each galaxy are 
of the same metallicity. Total mass-loss rates, determined 
from CO emission, may provide a better constraint of the 
gas density. This would require observations with the full 
ALMA array for objects at the distance of the Magellanic 
Clouds (Woods et al. 2012, MNRAS, in press), however, the 
correlation between CO and total gas mass is also uncertain. 



5.2 Onset of crystallinity 

Theore tical dust conden s ation studies pr edict th a t the 
dust dSogawa fc Kozasal Il999h or gas (|Tielens! Il998l ; 
iGail fc Sedlmavrll 19991 ) column density in the circumstellar 
envelope is critical for the formation of crystalline silicates. 
Our observations show that, at high dust densities (Mdust > 
3 x 10~ 8 Mq yr _1 ), the majority of the sources in our sample 
exhibit silicate grains with a crystalline component. How- 
ever, we first detect evidence for crystallinity at a low dust 



1 M 



©yr 



column density corresponding to Mdust ~ 5 x 10" 
Approximately 10 per cent of the sources with dust mass- 
loss rates between 10 -9 ' 3 and lO _8 M0yr _1 show the dis- 
tinctive structure of crystalline silicate bands. A transition 
in the silicate structure is seen when dust mass-loss rates of 
10 -8 to lO _7 M0yr _1 are reached. Here, the silicate grains 
undergo a transition from only a few per cent of sources 
displaying crystallinity to over 80 per cent of the sources, 
in less than 1 dex in dust mass- loss rates. This transition 
appears to be sharper for the dust mass-loss rates compared 
to the gas mass-loss rates, although the uncertainties are 
greater for the gas mass-loss rates. Although tentative, this 
suggests that crystallinity is predominantly correlated to the 
dust mass-loss rate, and thus annealing of amorphous sili- 
cate grains by radiation is probably the primary formation 
mechanism for crystalline silicates in the outflows of AGB 

and RSG stars. 

Alternatively, ISloan et all (|2008l ) argue that in low- 
density winds grain growth occurs over a long period, 
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allowing the atoms to arrange themselves in the most 
energetically-favourable lattice structure before being locked 
in place as the grain accumulates more layers. At higher den- 
sities grains will accrete material faster, preventing the mi- 
gration of atoms. If this was the case one would most likely 
expect to see a bimodal distribution of crystalline sources, 
which is not supported by our results. 

For sources undergoing the most extreme mass-loss 
rates a combination of crystallisation processes may be in 
effect. Here conditions are favourable for both thermal an- 
nealing of the amorphous grains and crystalline grains form- 
ing by direct condensation from the gas phase, where the 
higher densities in the dust-formation zone raise the sili- 
cate condensation temperature above the glass t emperature 
l|Lodders fc Feglevll 19991 ; iGail fc Sedlmavrll 19991 ). This com- 
bined effect may in part explain why the sources with high 
mass-loss rates tend to have stronger crystalline features. 



5.3 Mineralogy and metallicity 

Our results provide strong evidence for a change in crys- 
talline silicate production with metallicity; here we explore 
some of the possible reasons for these differences. 

The molecules available for dust production in O-rich 
stars are limited by the abundances of heavy elements which 
are not produced on the AGB. The abundance of Mg, Si, Ca, 
Al, and Fe in the photosphere reflect the initial abundances 
of the molecular cloud from which the star formed, while the 
surface abundances of C and O are altered by dredge-up and 
hot bottom burning. Dredge- up only increases the [O/H] 
abundance by small amounts. However efficient hot bot- 
tom burning in low-metalli city, high-mass AGB stars may 
result in [O/H] depletion l|Ventura et alj [2002). For stars 
with solar-like abundances, both Mg and Fe are more abun- 
dant than Ca and Al (Mg/Fe being approximately unity). 
These abundances change with overall metallicity. For non- 
q-capture products like Al, this effect is even more severe 
l|Wheeler et al.lll989l ). 

For solar-metallicity stars, the classical dust condensa- 
tion sequence in the outflows of O-rich evolved stars pre- 
dicts that aluminium- or calcium-rich dust grains form first, 
at temperatures around 1400 K, before silicon and magne- 
sium condense i nto silicates at slightly l o wer temperatures 
llTielensI ll99Ct lLodders fc Feglevl Il999l ; IGail fc Sedlmavrl 
1999). A parallel condensation sequence, involving mag- 
nesium and silicon, starts with the direct condensation of 
forsterite from the gas phase at temperatures around 1050 
K. This is subsequently transformed into enstatite through 
reactions between forsterite and Si02 gas, according to the 
equilibrium reaction: 



2MgSi0 3 Mg 2 Si0 4 + SiQ 2 . 



(5) 



Dust condensation for equilibrium conditions has been 
suggested to follow the sequence below: 



AI2O3, Alumina 
Ca 2 Al 2 Si07, Gehlenite 

a- 

CaMgSi 2 6 , Diopside Mg 2 Si0 4 , Forsterite 

MgAl 2 O4 , Spinel MgSi0 3 , Enstatite 

CaAl 2 Si 2 8 , Anorthite (Mg, Fe) 2 Si0 4 , Olivine 

The condensation sequence in O-A GB stars is s l ightly 
different from that of supergiants (|Speck et al ] l2000l : 
IVerhoelst et al.ll2009h . however, the end-point for both se- 
quences results in the formation of magnesium-rich silicates. 
The lower detection rate of crystalline features in RSGs may 
be a consequence of this alternate condensation pathway, 
due to the delayed production of Mg-rich amorphous sili- 
cates. We speculate that this production is a necessary pre- 
cursor for the formation of forsterite and enstatite, as it 
provides suitable material from which to anneal crystalline 
silicates. 

The crystalline silicates in our sample are comprised 
of Mg-rich olivines and pyroxenes (forsterite: Mg 2 Si04 and 
enstatite: MgSiOs). Recently, the crystalline features in the 
ISO spectra of RX Lac, T Cep, T Cet and R Hya were 
found to be bett er matched w i th th e iron-rich silicate: 
Mg i 8 Fei. 82 Si0 4 (jPitman et all l2010l : iGuha Nivogi et all 
2011). These sources occur early on in the silicate emis- 
sion sequence (which characterises the shape of the 10 /im 
complex; see ISloan fc Pried [l995) and lack the classical 
emission/absorption features which have been attributed to 
amorphous silicate dust, they also tend to have optically 
thinner dust shells. However, Herschel observations of the 
69 fim crystalline silicate band in more evolved AGB stars 
do not find eviden ce for iron in the lat tice structure of the 
crystalline olivine (|de Vries et al.ll201ll ). The sources in our 
sample lie later on the in the sequence, these tend to have op- 
tically thicker shells and a greater abundance of amorphous 
silicates. Consequently, Fe-rich crystalline silicates may only 
form early on in the AGB, in low-density winds. 

At solar metallicities, forsterite is the dominant crys- 
talline species, but as the metallicity drops enstatite be- 
comes preferentially formed. We consider five possibilities 
for this difference below. 



(i) The elemental abundance of Mg, Si and O during dust 
condensation will also influence the forsterite-to-enstatite 
ratio. If either Mg or O are depleted then enstatite could 
preferentially form. One might expect the abundances of Si 
and Mg to scale from the solar values by the same amount, 
in which case a lower relative abundance of O would be 
required for enstatite production. This could be due to wa- 
ter formation locking up free oxygen, changes in the natal 
C/O ratio with metallicity, or efficient hot bottom burning 
in low-metallicity high-mass AGB stars. 

In more metal-poor galaxies, dredge-up of newly produced 
carbon has a great er effect on the C/O ratio in the photo- 
sphere of the star (|Lattanzid [2003). A higher abundance of 
CO reduces the the amount of free oxygen available for dust 
production, which might result in the preferential formation 
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of grains with a pyroxene stoichiometry as they require less 
oxygen. 

The chemical abundances in massive AGB stars can also 
be modified by hot bottom burning. If the temperature's are 
high enough (T > 3 x 10 6 K) for the NO cycle and Ne-Na 
cycle to become active, oxyg en will become dep leted and 
sodium should be enhanced (|Caciolli et al.l l201ll ). Oxygen 
depletion (and associated sodium enhancement) has been 
observed in several low metallicity g lobular cluster stars 
|Kraft et al.ll 19971 ; [Cratton et al.ll200lP ). These stars have in- 
sufficient mass to have undergone hot bottom burning, how- 
ever, it may be possible that the abundances of the primor- 
dial gas was enriched in He, C, N and Na, but depleted in O 
by an early population of m assive AGB stars (jVentura et al.l 
l200ll ; lDecressin et alj|2007h . 

(ii) Differences in the temperature gradient in the out- 
flows of O-AGB and RSG stars in the Galaxy, Magellanic 
Clouds and Galactic globular clusters, will alter the position 
of equilibrium in the forsterite/enstatite reaction. As the 
temperature of the circumstellar envelope decreases, the po- 
sition of equilibrium in Eq. 0will shift to favour the exother- 
mic reaction, in this case the production of enstatite. Con- 
sequently the MgSi03/Mg2Si04 ratio would increase. This 
reaction requires that equilibrium is reached. 

(iii) Alternatively, non-equilibrium conditions during 
dust condensation may alter the order that olivines and py- 
roxenes form. In denser environments, the dust shells are 
more likely to reach equilibrium due to the increased prob- 
ability of interactions with other atoms/molecules, thus for 
Galactic sources where the partial pressure of dust is great- 
est, the circumstellar envelopes are more likely to reach equi- 
librium than sources at lower metallicities. Under equilib- 
rium conditions it is predicted that forsterite will form before 
enstatite. If forsterite condenses first then Mg will deplete 
at twice the rate of Si. Conversely, if enstatite condenses 
first the [Mg/Si] ratio will remain constant until the Si is 
exhausted. 

(iv) Furthermore, the higher partial pressure of dust- 
forming elements in the envelopes of Milky Way sources 
likely results in the growth of larger grains. The size of 
forsterite grains may also be an important factor in enstatite 
production: for large grains the surface area per unit volume 
decreases, reducing available reaction sites for the infusion 
of Si02 in the crystalline lattice. This curtails the amount 
of forsterite available for conversion to enstatite, and may 
result in large forsterite grains coated by a enstatite mantle. 

(v) Nucleation of Mg, SiO and H2O in the conditions 
relative to stellar outflows could also plausibly form en- 
statite l|Goumans fe Bromlevll2012T ). For silicates to nucle- 
ate from gas phase molecules, they must follow a series 
of thermodynamically- favourable reaction pathways, re- 
sulting in a small metastable silicate cluster with an en- 
statite stoichiometries. Once a cluster has formed these must 
grow under appropria te temperature cond itions to become 
macroscopic particles |Gauger et al.|[l99(t ). The subsequent 
growth of this cluster is extremely exothermic, this increase 
in internal grain temperature may be sufficient for the par- 
tial crystallisation of the forming dust particle. Why this 
nucleation pathway would be favoured in the metal-poor 
environments of the Magellanic Clouds remains unclear. 



6 CONCLUSION 

We have analysed Spitzer and ISO infrared spectra of 217 
O-AGB and 98 RSG stars in the Milky Way, Magellanic 
Clouds and Galactic globular clusters, to explore the onset 
of crystallinity and investigate how the mineralogy depends 
on the physical and chemical conditions of the star's enve- 
lope. Dust mass-loss rates were established through spectral 
energy distribution fitting with the GRAMS model grid and 
the mineralogy of the crystalline features determined from 
the spectra. The main results of this study are summarised 
as follows: 

• We detect crystalline silicates over 3 dex in dust mass- 
loss rate down rates of ~10 -9 Mq yr _1 . 

• Crystalline silicates are more prevalent in higher mass- 
loss rates objects, though sources undergoing the highest 
mass loss do not show the 23-fim forsterite feature. This 
is due to the poor contrast of the low temperature of the 
forsterite grains above the continuum and in some cases 
(self-)absorption of the short wavelength (A < 25/rni) crys- 
talline silicate features. 

• The dust mass-loss rate appears to have a greater in- 
fluence on the crystalline fraction than the gas mass-loss 
rate. This may suggest that the annealing of amorphous sil- 
icate grains by radiation is probably the primary formation 
mechanism for crystalline silicates in the outflows of AGB 
and RSG stars. 

• O-AGB stars have a higher proportion of sources with 
crystalline silicates features than RSGs, however, there is 
little variation in the structure of the crystalline silicate dust 
for O-AGB and RSG stars. 

• We report the presence of a newly detected 22.8 /im 
emission feature in the spectra of Milky Way AGB and RSG 
stars. 

• We detect a change in the crystalline silicate mineral- 
ogy with metallicity, with enstatite seen increasingly at low 
metallicity, while forsterite becomes depleted. This variation 
in the olivine-to-pyroxene ratio can be explained by a num- 
ber of possible mechanisms. 
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